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Abstract 

We investigate baryogenesis at a first order electroweak phase transition in the presence 
of a CP violating condensate on the bubble walls, in the regime in which the bubble walls 
are 'thick', in the sense that fermions interact with the plasma many times as the bubble 
wall passes. Such a condensate is present in multi-Higgs extensions of the standard model 
and may be formed via an instability in the minimal standard model. We concentrate 
on particles with typical thermal energies in the plasma, whose interactions with the wall 
are accurately described by the WKB approximation, in which a classical chiral force is 
evident. The deviations from thermal equilibrium produced by motion of the wall are then 
treated using a classical Boltzmann equation which we solve in a fluid approximation. 
From the resulting equations we find two effects important for baryogenesis: (i) a classical 
chiral force term due to the CP violating background, and (ii) a term arising from hyper- 
charge violating interactions which are pushed out of equilibrium by the backgroundfield. 
Provided the wall propagates slower than the speed of sound, both terms lead to the dif- 
fusion of a chiral asymmetry in front of the wall. This can produce a baryon asymmetry 
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of the observed magnitude for typical wall velocities and thicknesses. 
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1. Introduction 

In this paper we present a detailed discussion of electroweak baryogenesis induced by a 
CP violating condensate field on 'thick' bubble walls during a first order electroweak phase 
transition. 'Thick' in this context, and as we will see more precisely, means that the mean 
free time for a fermion propagating in the plasma is short compared to the time taken for 
the wall to pass. In this case one expects that the non-local quantum reflection effects, 
which such a CP violating condensate has been previously shown to produce [1],[2],[3], 
will be suppressed due to scattering. Instead we look for local classical effects which can 
produce a driving force for baryogenesis. A classical treatment has, as we shall see, many 
advantages in that there is a systematic framework (a Boltzmann transport equation) 
within which to compute the nonequilibrium effects in which we are interested. A shorter 
version of this work has already appeared [4]. 

There is still considerable uncertainty as to what the relevant bubble wall thickness 
and speed are. Calculations are difficult [5] and strongly dependent upon the still poorly 
determined effective potential. A recent detailed study by one of us (T.P.) and G. Moore [5] 
using some of the methods developed in this paper indicates, within the minimal standard 
model, for Higgs masses of order 30-70 GeV, and ignoring possible nonperturbative effects, 
a wall velocity ~ 0.4 and a wall thickness of order ~ 25T -1 . Thus typically quarks interact 
very many times via gluon exchange processes as they cross the wall. If this is indeed the 
relevant regime, then for top quarks at least (the most obvious mediator of electroweak 
baryogenesis since they couple most strongly to the bubble wall if one has standard model- 
like Yukawa couplings) , the particle- wall problem cannot be treated without including the 
effects of strong (QCD) scattering from the plasma. 

In an accompanying paper [6] we have introduced the essential ideas motivating the 
calculations which we undertake here. We showed there that the Lagrangian for a fermion 
propagating in the background of a bubble wall in a two Higgs doublet extension of the 
standard model can be written 

C = * 7 ^(<9 M - z<M^ 7 5 )* " mm (1) 

where g A Z^ = g A Z^ T - l^rp^d/A g A = ±\g, g = a/ ' g\ + g\ [6]. The + sign is for 

up-type quarks and (left-handed) neutrinos, the — sign for down-type quarks and charged 
leptons. g\ and gi are the gauge couplings of the SU(2) and U(l) gauge fields; v\ and 
v 2 the magnitudes of the vevs of the two Higgs doublets, the first of which is taken to 
couple to the fermions through Yukawa terms. The two contributions to Z M come from the 
CP-odd scalar field 9 which is the relative phase of the two Higgs fields, l -p\ l -p\ = Re 10 , and 
the condensate discussed in [7] which may be present even in the minimal theory. All 
the vector couplings to Z are removed by using the remaining unbroken vector symmetries 
to remove the pure gauge Z (the Z condensate piece is also pure gauge if we treat the wall 
as planar, and assume it has reached a stationary state in which the Higgs and gauge fields 
are functions of z — v w t). When the vevs vanish the remaining axial Z is a pure gauge field 
(for both fermions and Higgs fields, with charges g A = and g A = —\g for the charged 
and neutral Higgs components respectively, since g A = [|(?3 — Y) + \ {B — L)] g [6] ). 

However, on the bubble wall this pure gauge field has as we shall see very tangible 
CP violating effects, even on particles with typical thermal energies. The axial gauge field 
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condensate Z formed on the bubble walls at the first order phase transition violates CP 
spontaneously. Z° is odd under CP, and the spatial components Z % are CP even so that a 
bubble on which the spatial vector Z 1 points out everywhere is mapped under CP to one 
on which it points in. (Thus in the latter case it is actually the gradient Z' which violates 
CP). 

The departure from thermal equilibrium which this brings about as the bubble wall 
moves through the plasma will source baryogenesis. In the companion paper [6] we dis- 
cussed the case originally investigated by Cohen, Kaplan and Nelson (CKN) [1], when 
the fermion is treated as free on the bubble wall. The conditions for the validity of this 
treatment are discussed in [6], and roughly require that the wall thickness L is much less 
than the mean free path of the fermion A/. In this paper we consider the perturbations in 
the plasma produced by the CP violating background in the possibly more realistic regime 
of wall thickness for which the fermions interact frequently on the bubble wall, a condition 
which will be specified more precisely in the course of our treatment. 

As we have discussed in [6] one might expect that the inclusion of interactions would 
wipe out any interesting CP violation if the effect is a non-local quantum mechanical one. 
This is precisely how the thin wall limit has been understood [1], [2], [3]. However as we 
noted in [6] the WKB limit is not as trivial as it appeared when viewed simply in terms 
of reflection coefficients for monotonic wall ansatzes. The fact that the dynamics of WKB 
particles are non-trivial - in particular that WKB particles propagate like particles in a 
classical CP violating potential - suggests that there may be interesting effects which do 
survive when the scattering on the wall is included. 

Secondly, as noted originally by CKN [8], such a CP violating background perturbs 
the energy levels of particles and can push processes out of equilibrium locally. The 
original form of this 'spontaneous' baryogenesis took this perturbation to the energy to 
be modeled by a fermionic hypercharge potential, and calculated the resultant chemical 
potential driving baryon number violation by imposing constraints on exactly conserved 
quantum numbers. Both of these aspects of the calculation have been criticized. Dine and 
Thomas [9] pointed out that the fermionic hypercharge potential cannot be appropriate 
as the effect does not vanish as the Higgs vev vanishes. Subsequently, we pointed out [10] 
that the imposition of the constraints neglects transport processes which tend to restore 
the region to a local thermal equilibrium in which there is no baryon number violation. 
The treatment we will present here will take account of both these criticisms and show 
that the essential effect does survive and can also, when transport is taken into account, 
source significant perturbations in front of the bubble wall [11]. This has also been pointed 
out in a recent paper by CKN [12]. One of the objectives of the present work is to bring 
these previously unconnected pieces into one coherent framework which includes all the 
important effects and clarifies their relation to one another. 

This paper is organized as follows. Section 2 discusses the WKB treatment of parti- 
cle dynamics described by the Lagrangian (1). In section 3 we introduce the Boltzmann 
equation and discuss the fluid approximation with which we truncate it to an analytically 
tractable form. In the following section we derive the resulting fluid equations. In section 
5 we analyse these equations in several steps, illustrating how perturbations may be gen- 
erated in front of the wall and identifying the parameters which determine the behavior 
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of the solutions. We derive the reduced equations needed for the calculation of baryon 
production in much of the favored parameter space of wall thicknesses and velocities. In 
section 6 we analyze these equations, treating the two source terms separately, and calcu- 
late the resulting baryon asymmetry in each case. In section 7 we compare the results we 
have found to those obtained in the thin wall case in [6]. In section 8 we conclude with a 
summary of the paper and a discussion of its shortcomings and directions for future work. 

2. WKB Dynamics 

The WKB approximation to the dynamics of particles in the background of the bubble 
wall is good provided the length scale on which this background varies is long in comparison 
to the de Broglie wavelength of the typical thermal particles we wish to describe. This is 
simply the requirement that the thickness of the bubble walls L be greater than T -1 . As 
we have indicated above, this is a very reasonable expectation. 

To describe the WKB 'particles' we turn to the Dirac equation derived from (1). The 
dispersion relation is obtained as follows. In the rest frame of the bubble wall we assume 
that the field Z^ = (0, 0, 0, Z(z)), and we can boost to a frame in which the momentum 
perpendicular to z is zero. In this frame the Dirac equation reads (after multiplying 
through by 7 ) 

ic>oV = 7 (-H 3 ^ + m)^-^ZE 3 V (2) 
where E 3 = 7°7 3 is the spin operator. Setting ip ~ e -tEt+i Pz z^ we gee ^ e energy is 
given by the usual expression for a massive fermion plus a spin dependent correction. The 
eigen-spinors are just the usual free Dirac spinors. Returning to the p± 7^ frame amounts 
to replacing E with \/E 2 — p 2 ±: from which we find the general dispersion relation in the 

wall frame, 

E=[pl + (^pJT^T9AZ) 2 ]^ E 3 = ±l (3) 

where E 3 is proportional to the spin S z as measured in the frame where p± vanishes. The 
same dispersion relation holds for antiparticles. The particles we are most interested in for 
baryogenesis are left handed particles (e.g. ti) and right handed antiparticles since 
these couple to the chiral anomaly. Note that they couple oppositely to the Z field. 

In Figure 1 these dispersion relations are plotted for p± = for (i) m > g A Z ', (ii) 
m < gAZ and (iii) m = 0. We see how the branches deform into one another as we turn 
on the mass. In particular we note how the left- (L) and right-handed (R) branches break 
up into two pieces and form the E 3 = ±1 branches, as the mass couples the two chiral 
components on the wall. Correspondingly it is straightforward to see how the eigen-spinors 
(the usual Dirac massive eigen-spinors) become chirality (7 5 ) eigenstates for p z » m. The 
conservation of spin on the wall gives a simple picture of how this happens - an in-going 
left-handed particle incident on the wall evolves in an eigenstate of spin on the wall. If its 
momentum is reversed it emerges as a right handed particle since its spin is conserved. 

For what comes below we will find it useful to keep this identification in mind. We 
label the states on the wall by the states they deform into as the mass is turned off and 
write the dispersion relation as 

E L ' R =\pl + (VpI + m 2 ± sign(p z )g A Z) 2 ] * 

1 ( 4 ) 
[Pi + (Pz ± 9aZ) 2 ] 2 as m -> 
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where the +(— ) singles out the states which become L (R) in the unbroken phase. (We 
assume the wall propagates from left to right, so that incident particles from the unbroken 
phase initially have p z < 0.) The anti-particle L of a left-handed particle L is right-handed 
and so we write the dispersion relation for anti-particles as 

E L ~ R = vb R = [Pi + (VpJ+^T sign(p z )g A Z)^ . (5) 

In the WKB approximation we take each particle to be a wavepacket labelled by canonical 
energy and momentum (E = p°,p) and position x. To compute the trajectory of such a 
wavepacket, we first calculate the group velocity 

Vi =Xi = d Pt E (6) 

and second, using conservation of energy E = i-idiE + pid Pi E = 0, we find 

Pi = -diE (7) 

Together these constitute Hamiltons equations for the particle. The momentum of the 
particle is not a gauge invariant quantity, but the particle worldline certainly is, and we 
can for example calculate the acceleration from the Hamilton equations. We find 

dv z 1 d z (m 2 ) d z (qAZm 2 ) , . 

where of course E and p± are constants of motion. The first term describes the effect of 
the force due to the particle mass turning on, the second the chiral force. In the massless 
limit the latter vanishes, as it should because in this case the chiral gauge field can just be 
gauged away. 

We shall later need a few more explicit expressions, namely 



P± 1 f , ry \Pz\ 1 

v ± = ~^ v z = -[p z ±g A Z-===\ 



E E y/pl + 



m 2 " 



. y/pl +m 2 ± sign(p z )g A Z f . d 2 
P±=0 Pz = = [±sign{p z )gAO z Z + 



(9) 



m 



E — 2^pJ+ 



in 



2 - 



As m — > we recover the equations of motion for the L and R particles in the pure gauge 
field Z. 

Before we move on let us remark on our neglect of one effect in deriving these dispersion 
relations - thermal corrections to the fermion self-energies. In Appendix C we discuss how 
the dispersion relations and consequently the force terms are modified when one includes 
this effect. Depending on whether the Z field is 'strong' or 'weak' in the sense defined in 
the Appendix, one can show that the resulting dispersion relations are either those given 
here plus small corrections, or slightly modified ones which lead to force terms with altered 
momentum dependence. This is a complication we will ignore in the rest of this paper. 
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3. The Boltzmann Equation in the Fluid Limit 

Since most particles in the plasma are well described by the WKB approximation with 
respect to their interactions with the wall, we proceed to treat the fluid of such excitations 
as consisting of classical particles with definite canonical position and momenta, and energy 
given by the derived dispersion relations. The Boltzmann equation for the phase space 
density f(p, x, t) is 

d t f = d t f + Z-d s f + p-d p f = -C[f]. (10) 

The collision integral C[f] describes how the phase space densities are changed by interac- 
tions. The dominant interactions which we will consider are, because of Debye screening, 
short ranged, and to a first approximation may be treated as pointlike. They are to be 
calculated at a given spatial point x using the Dirac spinors appropriate to the local value 
of the background fields m(x), and Z(x), taken to be constant. The Boltzmann equation 
is in principle solvable, but in order to make it analytically tractable we shall we consider 
an approximation (truncation) which we expect to be quantitatively reasonable based on 
a perfect fluid form for the phase space density. If interaction rates are fast, the collision 
integral forces the phase space density towards a form which minimizes it, namely local 
thermal equilibrium. For the case of a single fluid this means the form 

f& X ' = e f3(~f(E-v.p))-v) ± I ( n ) 

where /?, p and v are functions of x and £, and 7 = 1/(1 — -u 2 ) 1 / 2 . Imposing this form 
on the left hand side of the Boltzmann equation and integrating to find the moments, 
one arrives at a set of coupled equations for the functions (3(x, t), p(x, t), v(x, t) which 
parametrize the phase space density. Expressed in terms of the energy density p, number 
density perturbation n, pressure p and velocity v these become the familiar fluid equations 
in the case of a free fluid: 

d t n + V • [ynv\ = 

dtp + V- [ip + phv] =0 (12) 

d t [(p + phv)]+Vp = 0. 

We are going to treat the approximately left-handed excitations L and their antiparticles 
L as two fluids, making an ansatz of the form (11) for each. 

This requires some justification, because the dominant interactions which bring our 
fluids to the form (11) are the same ones which damp away the temperature and velocity 
differences with respect to the background. This is not true of the chemical potentials 
which are only attenuated by slower chirality changing processes. We keep the temperature 
fluctuation 5T in order to quantify the rate at which the form (11) is approached. We use 
the velocity perturbation in the form (11) to model the anisotropic response to the force. 
We will in fact see that the precise form of this perturbation does not enter the final result, 
and its only role is to allow the particles to move in response to the force, and set up a 
chemical potential perturbation. 

In the present case we wish to determine how different species are perturbed, by source 
terms which will enter these equations in a way we will calculate. In particular we must 
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distinguish between particles and anti-particles as it is the difference in the perturbations 
to these which is needed to source baryon production through sphaleron processes. To do 
this we will treat each particle species as a fluid described by the fluid ansatz (11), (in 
which the functions /?, fx, v are allowed to be different for each species) which makes the 
self-interaction collision terms zero. We then superimpose on this the interactions of these 
fluids, which will lead to terms in the equations damping all these perturbations to the 
local thermal equilibrium for the whole fluid. 

Is this a good approximation? There are three conditions which must be fulfilled: 

(i) The interactions must be fast enough to ensure the system is maintained in the 
approximate form (11) as the wall moves. This should be a good approximation if the time 
the wall takes to pass is long in comparison to the time scale r for the system to attenuate 
fluctuations away from this form i.e. we want — > r. What should we take this time 
scale to be? We will see in due course that this question receives an answer within our 
calculation: We should take the scale r to be approximately D, the diffusion length. 

(ii) The rate at which the system is brought to this local equilibrium form should 
be faster than the rate at which the perturbations we keep in our ansatz are attenuated. 
This is not the case, because the same gluon exchange processes which force / to the form 
(11) also damp away the velocity and temperature perturbations. Thus (11) is unlikely to 
be a very accurate description of the precise form of those perturbations. The chemical 
potential fluctuation will be the crucial term in our answer, which shall determine the final 
baryon asymmetry, and this is attenuated only by processes which change the number 
of particles in our 'fluids'. As we will discuss at length these are indeed typically much 
slower than the gauge boson exchange processes which damp ST and v. The temperature 
perturbation shall actually play little role - we keep it merely to see how thermalization 
occurs and to estimate its rate. 

(iii) The mean free time for particle interactions should be long compared to their 
inverse energies, in order that we can descibe the physics in terms of a set of 'free particle' 
eigenstates of well defined energy. This condition means the 'width' of a state AE should 
be much smaller than its energy E. The particles of central interest in this paper have 
energies E > T, and AE ~ g 2 T, so the condition is reasonably well satisfied even when 
one includes the strong interactions. 

We should also mention a subtlety about the distribution (11). The form is dictated 
by the quantities which are conserved in the local interactions. That these are E = p a 
and p is a non-trivial fact in this case in which the gauge symmetry is broken. In the 
presence of a pure gauge field, it is not immediately obvious whether one should take just 
the canonical momentum, or perhaps some linear combination with the gauge field, as the 
correct conserved quantity. 

For a conserved gauge charge, which one takes is immaterial since the difference can 
always be absorbed in the chemical potential, which is determined by the condition on the 
conserved charge. But here the gauge symmetry has been broken, and we thus need to 
determine what the correct conserved energy and momentum are in all local interactions, 
including those which violate the broken gauge charge, which we assume to take place in 
approximately constant background fields. By a simple calculation one can show that the 
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'canonical' stress tensor 



0^ = ^id v i) 



(13) 



is in fact conserved in constant background fields, since from the Dirac equation derived 
from (1), 

d^ v = {d v m)W + <7 A (d M ^)V^ 7 5 V> (14) 

The conserved energy and momentum are then easily calculated from the eigenstates we 
discussed in the previous section, and are indeed just the canonical energy p° = E, and 
momentum p. This is as to be expected from translation invariance. 

We also need to specify what precisely the fluids are as this is ambiguous once we 
turn on the background in which the dispersion relations and particle eigenstates become 
different. In the unbroken phase we take the fluids to be the chiral eigenstates. Which 
states described by the dispersion relations in Figure 1 do we take to make up the fluids 
described by the phase space density in (11)? The answer is that foreshadowed by the 
way we wrote the dispersion relations in (4). The motivation for this choice is that most 
particles in these fluids are then to a good approximation in a given chiral state. This of 
course breaks down completely for low momentum states p z ~ m which are effectively equal 
mixtures of the chirality eigenstates. By dividing the fluids in this way we also misdescribe 
the dynamics of these low momentum states by connecting the wrong branches at p z = 0. 
Most particles in the fluid will not change direction between scattering and so this should 
be a small effect. We will return to these points later. The essential point is that the 
effects we describe will be dominated by particles at thermal energies. The force is felt by 
all particles in the plasma, not just those at low momentum. 

4. Fluid Equations 

We now proceed to derive in detail the truncation of the Boltzmann equations which 
we have just described. The substitution of the ansatz (11) in the left hand side of (10) 
gives, in the rest frame of the wall in which the energy is time independent, 

- [pdtH+Pp z dtv-dtP{E-n-vp z )\ /'- [pd x n+Pp z d z v-d z p{E-n-vp z )] f'd Pz E+[3vf'd z E 

(15) 

where /' = ^( e Jj-i ) and x = (3(E — ft — vp z ). We take v = v z and work to leading order 
in v. Making the substitution ft = fx ± vqaZ (with the sign chosen appropriately for the 
left-handed and right-handed particle and antiparticle fluids ) and writing p z ±gAZ = k z , 
the physical momentum in the zero mass limit, we then have 

- [Pd t fi + [3k z d t v - d t p(E -ft- vk z )} f - [(3d z ji+(3k z d z v - d z (3{E -ft- vk z )} f'd Pz E 

+(3vf'[d z E T 9Ad z Z.d Pz E}. 

(16) 

We make this change in variables so that our equations will be in the familiar form (12) 
as m — > 0. The variable ft is really what one would usually call the chemical potential 
in the m = limit, as fluctuations in number density are proportional to it. The fi in 
(11) is in fact a gauge dependent object and ft is simply the appropriate gauge invariant 
chemical potential which appears in real physical quantities calculated from the distribution 
function. 
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The crucial term here is the coefficient of the velocity v. We are working in the rest 
frame of the wall and so write v = —v w + v since we are interested in the perturbations to 
the background which is a plasma moving by the wall at velocity —v w in this frame. Then 
we see that there is a source term proportional to v w with coefficient 



p M^? ±n S n(p,) g ± Z [±siMa8mi _ \PA } + _^_] (17) 



-m 



where ± correspond to the signs in the dispersion relations (4) and (5) for the particles 
and anti-particles in our (approximately) chiral fluids. This source manifestly vanishes as 
m — > and when integrated will give us the force on the plasma when the wall is pushed 
through it. Before doing this integration of these equations we consider also the right hand 
side of the Boltzmann equation. 

The collision integral C[f] in (10) is equal to the rate of change of the phase space 
density f(p, x) due to collisions. Considering only processes in which there are two incom- 
ing (labelled as 1 and 2) and two outgoing particles (labeled as V and 2') we can write it 
in the following form 

processes^ 1 VwPz >J Pi> .P2 ,P 2 ' V 7 (18) 

V[fi] = [/i/ 2 (l =F =F fv) - fvMl =F /i)(l =F / 2 )] 

where the precise spinors used as eigenstates will determine the various normalizations in 
the integral, \M\ 2 is the matrix element for the process, and =F is for fermions (bosons). 
Taking the fluid ansatz (11) for each particle species and doing a perturbative expansion 
in (3 = (about j3 Q = ^-), ^- and v we get, to leading order, 

E 2(F \, n J \M\^)H\^ p )[hf 2 {l T fv){lTh')] 

(iyj 



x p 



J- n 



where fi etc. are now the unperturbed distribution functions at temperature T Q . We have 
performed this expansion in the plasma frame variables for reasons we will explain below 
when we come to integrate this expression. The normalization factor is just the energy in 
the wall frame. XX/^j) means a sum over chemical potentials with a positive (negative) sign 
for in-going (outgoing) states. In the temperature and velocity terms we have assumed that 
the in-going and out-going 1 and V are in the same fluid and the same of 2 and 2', since 
this is the case for the dominant scattering processes. For each term in brackets we must 
take the fastest process which force these fluids to the same thermal equilibrium. The ones 
that attenuate the temperature and velocity perturbations are gluon exchange diagrams 
(for quarks) or weak boson exchange (for leptons), shown in Figure 2. These processes 
do not contribute to the chemical potential damping however since they do not change 
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particle number. Examples of processes which do contribute are the helicity flipping gluon 
exchange process (as in Figure 2, but in the presence of a mass term) which can occur on 
the wall, and the Higgs mediated decay like that in Figure 3, as well as sphaleron processes 
(both strong and weak). 

We now integrate the Boltzmann equation over J ct 3 p, J ct 3 pE and f ct 3 pp z ( wall 
frame variables) to get three moments which give us three first order differential equations 
for the three functions %£-,v,-M- characterizing the perturbations in each fluid. In Appendix 

J- o J- o 

A the terms obtained by integrating the collision integral are analyzed. Considerable 
simplifications occur provided the integrations are done in the plasma frame variables 
in which f a has the standard (unboosted) form. Thus we take the linear combinations 

3 3 

J <f p(E + v w p z ) and J <f p (p z + v w E) of the latter two integrations which are integrations 
over the energy and momentum in the plasma frame (to leading order in v w ). The velocities 
which appears in (19) are those in the plasma frame, but the difference of velocities in this 
frame is equal to the difference of the velocities in the wall frame to leading order in v w . 
We work perturbatively in v, ■&- keeping the leading terms in these quantities and their 

J- o J- o 

derivatives, dropping next order corrections to these coefficients in ^ and ^. 
The equations which result are 

8T . ft . 1 , m 2 , - v-^ ft ft Z 

-Vwijr) -av w {—) +gf +gv w (jv) = -T^2^ — -T^2^(—±v w g A — ) 

-v w (—Y - bv w (jr) + W + hv w ( w ) = 

lo lo 6 _° (20) 

_rT E t~ ~ f !» S l" -1 ^ /Z^ ± v ^3A7fr) 

± o ± o ± o ± o 

,&T V , u Ji v -, , (g A Zm 2 )' ^ _ 

(— ) +b(—) -v w v ±cv w -3 = -T v ^v 

where a = |^-, b = c = 9 = W^, h = and ( n are the Riemann (^-functions 

((2 = -^p C3 = 1-202, (4 = 7r 4 /90). To derive these equations in this form we have expanded 
E, d z E, d Pz E around m = 0. Except the force terms all the terms on the LHS have just 
the free fluid coefficients - precisely what the equations (12) give if one expresses them in 
terms of 8T : fj,, v ( for m = i.e. p = |p). The positive sign in the force term applies to the 
L and R fluids, the negative sign to the R and L fluids; the mass is that of the appropriate 
fermion. 

The terms on the RHS require some explanation. We have written 

r,£^ + r M £(A±^4) 

J-o J-o J-o 

E f + K E f + r ', ± WA^r) (21) 

J-o J-o J- o J-o 

r -E^ 

where n a and p Q are the unperturbed number and energy densities respectively at T Q . 
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To arrive at the RHS of these equations we have used symmetry arguments to show 
that various terms are zero for tree level processes we are interested in. This is discussed 
in Appendices A and B. 

In the collision terms as written in (19) we see that the substitution of the variable 
v gives the same expression with v replaced by v as it is only the relative velocity which 
is damped by this term. However the substitution of p, = p =f v w9aZ ± vgAZ is very 
non-trivial. The latter piece gives only a higher order correction to our equations. But 
the term v w gAZ only drops out to give the same expression with \i replaced by p, if the 
charges qa on the external legs in the process sum to zero. In the limit that the vevs of 
the Higgs fields vanish Z is a pure gauge field for the fermions and Higgs particles, and so 
this cancellation occurs for all processes which conserve any linear combination of electric 
charge and hypercharge. Thus f M and T' are calculated from decay processes which con- 
serve hypercharge, T M and T'^ from processes which explicitly violate it, the latter picking 
up a net contribution proportional to v w Z. These latter rates, being hypercharge (and 

violating, are vev squared suppressed and therefore any baryon asymmetry produced 
by them will vanish as the vevs do. When the vevs are nonzero, the Z field perturbs hy- 
percharge violating processes out of equilibrium locally, while the hypercharge conserving 
processes remain in equilibrium, simply because they conserve the charge associated with 
the gauge field. These latter "see" Z as a pure gauge mode, not a real gauge invariant field 
which shifts the energies. This is precisely the type of effect which CKN called 'sponta- 
neous' baryogenesis. We will discuss in section 6.2 what light this treatment throws on the 
question raised by Dine and Thomas in [9] about how the background should be modeled. 

Before we move on to solve (20) for some specific cases we make a few general com- 
ments: 

(i) We have dropped all the time derivatives because we are interested in stationary 
solutions. If we wish to understand, for example, how the stationary solutions are set up at 
the time of nucleation, we include the terms which are obtained from the time derivatives 
in (16) after integration. 

(ii) We have dropped all long range fields which result from the perturbations. This 
amounts to neglecting the effect of screening of electric charge and hypercharge on the 
solutions we will study. 

These terms can easily be included (through an addition to p ) and give a term 
proportional to the field E in the third equation. Making use of Gauss' Law this can be 
written in terms of the perturbations, thus coupling the LHS of all the equations to one 
another (see [6] section 7) . We will comment further on this point in the conclusions. 

(iii) As v w — > the only solution unperturbed at infinity is the trivial solution f3 = f3 Q 
and v = p, = 0. Any perturbation which may source baryon number (or indeed push 
any process out of equilibrium) arises due to the motion of the wall. This does not mean 
however that this static equilibrium is the unperturbed one which pertains when the field 

t The fact that we refer to these processes as 'hypercharge violating' has no particular 
significance. They violate any linear combination of hypercharge and electric charge. In 
particular they violate the axial charge gA which gives the coupling of the CP violating 
condensate to the fermions and Higgs fields. 
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is not present. This is because the energy in the distribution function is the perturbed 
energy. For example if we calculate the number density J c? pf\n=v=o and subtract the 
true unperturbed number density we do not get zero if we integrate e.g. over states with 
p z > 0. This equilibrium has in it an excess of right-moving spin ^ particles and an equal 
and opposite underdensity of left-moving spin — | particles. The way to understand this is 
by analogy with an electromagnetic potential which is screened. The thermal equilibrium 
reached in the presence of such a potential has an overdensity of particles in proportion 
to their charge. Here the "screening" of the force induced by turning on this field has the 
effect of dragging in particles as described by this distribution function with v = fx = 0. 
This was the essential point made in [10], except that the energy perturbation was modeled 
(inappropriately) by a purely fermionic hypercharge potential. 

(iv) As the vevs of the Higgs fields vanish the only solution is again the trivial one 
(without time dependence) since both the mass and the rates of hypercharge violating 
processes go to zero. In this case the distribution functions do also describe the true 
unperturbed plasma since in this limit the dispersion relation approaches the pure gauge 
one in (4). 

(v) We will explain below that in the case that we neglect the temperature fluctuations 
the system is reduced to the first and third equation in (20). Dropping all the force terms 
and setting all the decay rates to zero, to order v w we obtain 

Kjr)" = -3a^r„(A)' (22) 

1 o 1 o 

which describes pure diffusion h = —DV 2 n. We can then read off the relation D = ^T" 1 . 

5. Solutions of Fluid Equations 

We now turn to the analysis of the fluid equations (20) , with the goal of understanding 
the baryogenesis which results from the perturbations they describe. As discussed in 
section 2 of the accompanying paper [6] the anomalous baryon number violating process 
is perturbed from equilibrium by a difference in the distribution functions of left-handed 
fermions and their (right-handed) antiparticles. As we wish to study baryon production 
we therefore take the difference of (20) for left-handed particles and their anti-particles, 
and get 

± o " J- o 

tt(^) - r; £(f ) - r ; E(^4) ( 23 ) 

J- o -Lo J- o 

We have compactified our notation: T M now includes all decay processes and T M denotes 
only the hypercharge violating ones. The force terms in the first two equations in (20) 
cancel out because the gradient in the real mass affects particles and anti-particles equally. 
The parameters ST, v, ft now represent the difference in these quantities for particles and 



,5T. . . a 1 . 

-v w {—)-av w {^-)'+-v' 

- v ^y- bv ^Y o y+r 

&' + b(^)'-v w v' + F(z) 



12 



antiparticles, and the force 



F = 2cv w g A 



(Zm 2 )' 



= An 



w 



{Zm 2 )' 



(24) 



We have used the fact that particles and anti-particles couple in exactly the same way in the 
gauge boson exchange diagrams to cancel out the temperature and velocity perturbations 
of the other fluids. This removes the sum in the Ft and T v terms. The counting factor 
which results - over particles and anti-particles of all flavors - has been absorbed in the 
definition of IV and T v . For quarks we show in Appendix A that 



We note that the relation D = b/3aT v from (22) then gives D ~ 5/T in very good 
agreement with the value calculated by a different method in [6]. The only coupling to 
perturbations in other fluids remains in the sums for the decay processes. There are two 
distinct sources in (23) for the perturbations: 

• The force terms F(z). 

• The hypercharge violating processes which are perturbed from equilibrium as the 
wall passes when v w Z ^ 0. 

As the equations are linear we can separate these sources and study them indepen- 
dently. 

We will not attempt to solve the full set of equations in complete generality for each 
of these two source terms. Firstly, we will limit our scope by considering only the case 
where the source terms directly affect the top quark. Both the classical force term and (we 
will see below) the spontaneous baryogenesis terms are proportional to the mass squared 
of the fermion for which (23) describe the particle minus anti-particle perturbations. We 
will thus work in the approximation that only the top quark Yukawa coupling is non-zero, 
which is good if the fermion Yukawa couplings are (as in the minimal standard model) 
proportional to their zero temperature masses. We will not consider the case emphasized 
in [6] , where in a two Higgs extension the tau lepton can have a Yukawa coupling as large 
as that of the top quark (and hence a comparable finite temperature tree level mass). 

We now discuss: 

(i) Dropping the temperature fluctuations, so reducing the equations to just two cou- 
pled equations for ft and v, which can in turn be written as two second order uncoupled 
equations for the these two variables. 

(ii) The solution of the reduced equations for the classical force source term with the 
simplification that we neglect all decay processes. 

(iii) The problem with decay processes included, for the classical force source term. 

(iv) The source terms from hypercharge violating processes and how they source baryo- 
genesis. 

5.1 Thermalization and Validity of the Fluid Ansatz 

In section 3 we explained that we do not expect the temperature fluctuations to 
accurately reflect the perturbations in the plasma, because the tree-level gauge boson 



r„ = 3r T ~ a 2 s \n—T ~ ^ 



(25) 
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exchange processes which damp them away are exactly the same processes which damp 
away the perturbations to the distribution functions which we neglected in taking the 
ansatz (11). We kept them in our ansatz however in order to have a quantitative measure 
of the validity of our approximation. It is this point which we first consider here. 
Neglecting all decay processes the equations (23) are 

ST' ft' 1 , 

-v w — ~ av w — + -v = 

-v w — -bv w — + -v =-r T — (26) 

^T' , W t W X 

— + &£- - v w v' + F(z) = -T v v . 

J-o J-o 

One can convert these to a pair of second order uncoupled equations for 8T and p, (by 
integrating the first equation directly to find v), and then solve them for various sources F. 
We will not present here the details of this calculation as the results are of no importance, 
except for the following point. Two parameters - XdL and XtL - enter in determining the 
behaviour of the solutions, where 

— "TT A T = (27) 

D b — a v w 

are the two roots of (26) and L is the thickness of the wall. The diffusion root Xd describes 
the diffusion tail in front of the wall, while Xt describes the decay of perturbations behind 
the source. X^L is simply the squared ratio of the wall thickness to the distance a particle 



with diffusion constant D diffuses as the wall passes (~ yDt ~ ^DL/v w ); this is the 
parameter which, as we discuss below, will characterize what we call 'good transport' 
or 'poor transport'. XtL is the ratio of the time of passage of the wall (L/v w ) to the 
mean free time for temperature attenuating processes (~ T^ 1 ) and it is a measure of how 
efficiently the temperature perturbations are damped on the wall. One finds that when 
XtL > 1 the temperature perturbations in the solutions to (26) are damped by at least 

relative to the chemical potential fluctuations. Examining (26) one can in fact see 
this damping directly in the second equation. Taking all the derivatives to go as we 
see immediately that the suppression follows. When one incorporates the decay processes 
a similar conclusion follows provided IV >> T M , which applies (see rates given below). 

Because this condition for the damping of temperature fluctuations is the same as that 
of the validity of our initial ansatz we always take it to apply and reduce our equations to 
the first and third in (23) with 5T set equal to zero. Noting the relation between Ft and 
T v from (25) and T v = j^D~ x from (22), the condition becomes 

L . . 

v w < — — thermalization. (28) 

Oi-/ 

Still neglecting the decay processes, the equations simplify further: 

Jl' 1 n 
-av w — + -v = 

, „ °(zm 2 y „ (29) 

b— - v w v + Av w — —3 — = -T v v. 

J-o J-o 
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where A = leg a- The solution to the first equation gives v = 3av w -M- (using the boundary 
condition that /i = v = far in front of the wall) and substituting in the second equation, 
again using the relation T v = -^D~ l from (22), we get (to leading order in v w ) 

D T +Vw T = -J F{z) (30) 

As discussed in section 3 v is, like ST, not in itself to be taken to accurately describe 
the perturbations in the fluid. This is the case because we would expect there to be other 
anisotropic (in momentum) components of a general distribution function which we have 
neglected in our ansatz, which will be damped away by the same (tree-level gluon exchange) 
processes at approximately the same rate as the perturbation parametrized by v which we 
have taken. We cannot however consistently set it to zero in (23); but as we now see v is 
indeed damped by v w relative to ft, so that it itself will not contribute (at leading order 
in v w ) to the biasing of the baryon number violating processes, and neither does it enter 
in (30) which determines ft. Its only role is to mediate the force to the chemical potential, 
and we assume that any other anisotropic component would have led to approximately the 
same result. 

The fluid equations (26) are calculated to leading order in v w . A fuller analysis 
incorporating all orders in v w can be performed and shows that the velocity at which the 
leading order analysis breaks down is the speed of sound v w ~ v s = 1/ v^3 in the plasma. 
If the wall moves faster than this there is no solution in front of the wall and perturbations 
cannot propagate into this region. We are interested in the case when perturbations can 
propagate in front of the wall where the anomalous electroweak processes are unsupressed. 
Thus we assume 

v w < v s = —= (31) 
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5.2 Classical Force sourced Perturbations without Decay 

We now solve (30) for a chosen ansatz for the force term in order to illustrate the 
behaviour of the perturbations they describe and to gain some simple intuition for the 
more complex case where we include the decay terms. We take the following 'ramp' ansatz 
for the source term: 

^„^! = (fo(,-f) -*<?<* (32) 



T% [0 otherwise 

where F a is the constant force on the wall, and as above A = 2cgA- (30) has then a 
particular solution which is non-zero only on the wall, and a homogeneous solution which 
is an exponential describing pure diffusion ~ e~~ z . We impose the boundary conditions 

B(f)|t: = D(A)|:£ = ™£ (33) 

1 Q 2 1 Q 2 

obtained by integrating (30) taking p, to be at most step discontinuous across the bound- 
aries (so that its integral is continuous). Solving, requiring ft to be finite at — oo, we 
get 

^-^\-£l + (£z + ^-* i ^ ) -*<*<* 04) 

O " I r _ V W T n /-, _% Tm ( 7 _i) T, 

y[e — L --^ T (l-e — L )]e ~ (z z>± 

These solutions for -M- are sketched in Figure 4 for the two cases (i) ^ « 1 and (ii) 
^fy 1 » 1, when the solutions on the wall can be written 



-TP T \ 1 v w L 1 (z~ 2") I -1 1 ( VujL_ \ 21 v w L , , q 



(35) 



D 



In the first case we see that the solution mimics the behaviour of the driving force. From 
(29) one can see that this behaviour is generic for << 1, for expanding in v w one finds 
at leading order the solution 

ft A Zm 2 . . 

v = o, Y = ' Vw ~b^rf (36) 

is non-zero only on the wall. What the parameter v w L/D tells us is how efficient the 
diffusion is in bringing us to (36). As discussed in section 5.1 it is simply a ratio of the 
thickness of the wall to the distance a particle diffuses as the wall passes. This is the 
parameter which defines 'good transport'. The corrections to (36) describe a diffusion tail 
in front of the wall which has amplitude 

P I _ _ 1 v w L F a L 

T Q k ~ 2 D b [6n 

The factor of two arises from the fact that the diffusion in front of the wall is driven by the 
average amplitude of the potential on the wall, since particles are in this regime "seeing" 
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the whole wall. This is confirmed (see examples later) by calculating with other ansatzes. 
We thus conclude that in more general we would find the solution in front of the wall 



A = + ^ „„£<^!> e -^-i) (38 ) 

where (..) means the average value on the wall. 

When ]h! fj L » 1, i.e. when the transport is 'poor' over the relevant timescales, the 
solution looks like that in Figure 4 (ii). The amplitude on the wall rapidly approaches that 
in front of the wall which is 

P ° F ° (39) 



T a 1 2 v w L b 

so that the integrated amplitude Jl p, is suppressed relative to the 'good transport' case 

by (^t;) 2 - It is this integrated amplitude of p, (recall that it is the difference of particle and 
anti-particle chemical potentials) in front of the wall which will be the effective driving force 
for the final baryon asymmetry, when we assume the baryon number violating processes are 
immediately switched off on the wall. We note that the solutions (34) have the property 
that 

f + OO 

jj(z)dz = (40) 

In fact this can be shown directly by integrating (30) once, taking the solutions to be zero 
at ±oo. The integral of the chemical potential is zero because there is no net particle 
minus antiparticle creation in the absence of the decay processes. Thus the integrated 
contribution in front of the wall exactly cancels that on the wall, the excess of particles 
pulled onto the wall in response to the force being exactly cancelled by a deficit in the 
diffusion tail in front of the wall. In terms of baryon production it follows that if the 
sphaleron rate were unsuppressed on the wall exact cancellation would occur between the 
production in front of and on the wall so the source would not bias net baryon production. 
The extent of the cancellation which can occur is a sensitive function of the wall profile, 
but unless Z condenses only where the vev is very small this will not be significant. We 
will always assume that the baryon number violating processes are only turned on in front 
of the wall. 

This property (40) allows us to read off simply from (36) the integrated amplitude in 
front of the wall when << 1 (to leading order in this parameter) as 



oo — 



Jx A f Zm z 

¥ = Vw J "W" (41) 

J o J wall 1 o 



where again A = leg a (c and b being the numbers defined after (20) ). This in fact also 
proves (38) since we know that the solution in front of the wall ~ e~~^ z . We note that 
the result is in this regime independent of D and L, and we will see that the factor of v w 
also drops out when we calculate the baryon asymmetry. We will defer doing so explicitly 
until after the next section in which we discuss incorporating the decay processes in this 
analysis. 
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5.3 Decay Processes and Baryon Production 

Putting back the decay rates our equations are 

. U . t 1 _. - v— v ft v— v Z 

-av w { — ) +-v = -T^2_^ — + T fl 2^v w g A — 

J-o <j J o J-o (42) 

6(|-) / + F(z) = -r t ,tJ 

where we follow again the conventions in (23). ft is again the difference in particle and anti- 
particle chemical potentials and the sums are over the external legs of the processes which 
change the number of particles minus anti-particles. We have a pair of such equations 
for each fermion of (approximate) chirality, where F = ±Av w (+ for the L — L 

perturbations, — for the R — R perturbations) and m is the mass of the fermion, which 
we have assumed only to be non-zero for the top quark. The equations are coupled only 
through the sums in the decay terms on the right hand side, and it is this coupling which 
we now consider. 

As before we can decouple the variables jl and v and concentrate on a set of second 
order equations for the Jjl: 

d T +Vw ¥ - T ^t = -J F ' {z) Vw9A Y Q (43) 

The decay rates here have absorbed a factor of 1/a relative to the rates defined in (21). 
There are a similar set of equations for the v. We will not consider these further since, as 
discussed in the section 5.1, v is attenuated relative to p, and hence also its effect on the 
baryon number violating processes which is what concerns us here. 

What are the decay processes? Recall that our particle and anti-particle states are 
those of our (approximately) chiral fluids. They are the WKB eigenstates of the dispersion 
relations (4) and (5), which become pure chirality states in the unbroken phase. 

In the broken phase they are helicity states which are mixtures of the two chiralities. 

It is thus natural that we divide the processes into (i) those which occur in both phases 
and (ii) those which occur only in the broken phase. 

We consider these in turn: 

(i) In the unbroken phase our states are exact chirality eigenstates. The only pertur- 
bative processes which change chirality are those involving Higgs particles, like that shown 
in Figure 3, in which a left-handed quark flips chirality when it scatters off a gluon and 
emits a Higgs particle. Since we are taking only the Higgs-top Yukawa coupling to be 
non-zero this induces a coupling only between left-handed top quarks, right-handed top 
quarks and Higgs particles. The gluon cancels out when we subtract particles from anti- 
particles because it is its own anti-particle. The only other processes changing the number 
of particles of a given chirality in the unbroken phase are the strong anomalous processes, 
and the weak anomalous processes which are responsible for the baryon production. The 
former couple left-handed quarks of all flavors to right-handed quarks of all flavors directly; 
the latter couple all the left-handed quarks to all the left-handed leptons. 
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(ii) In the broken phase all the processes in the unbroken phase (except the baryon 
number violating processes) are still present, but there are many additional processes 
which couple particles on a given branch of our dispersion relations to other particle states, 
because of the mixing of helicity and chirality states when the mass is non-zero. 

When we calculate interactions between particles in these states we find that they 
can couple to one another through a 'helicity-flipping' gauge boson exchange. This is also 
an example of what we called a 'hypercharge violating' process: if we identify the ingoing 
states by the hypercharge of the state they deform into in the unbroken phase, hypercharge 
is not conserved. We have evaluated the rate for this process in Appendix B. There are 
many other such processes e.g. involving Higgs fields. There are also many flavor changing 
processes mediated by W bosons, but these are zero in the approximation that only the 
top quark Yukawa coupling is non-zero. 

We now make three further simplifications. Firstly we will assume that the baryon 
number violating processes can be neglected, except in their role as the source of net baryon 
number. We will see below that the condition that this be true is (see also [6]) 

r s < § (44) 

where T s is the rate for electroweak anomalous processes. v^/D is just the rate of capture 
of a diffusing particle by the advancing wall (the inverse of the time it spends diffusing 
before capture). 

The second further simplification we make is to neglect completely the decay processes 
in the unbroken phase involving Higgs particles on the external legs. The reason we do this 
is found in the study of these processes in [6]- in the unbroken phase the equations (43) 
are precisely the same diffusion-decay equations obtained there. The change to the results 
when these decay processes were incorporated was found to be a minor numerical one, and 
we assume that the same will be true here. In short the reason is that these processes do 
not drive the quantity sourcing baryon number (i.e. left-handed fermion perturbations) 
to zero. In the limit where they are fast enough to equilibrate locally they simply lead to 
a redistribution of particles amongst the left-handed fermions, right-handed fermions and 
Higgs particles. 

The final simplifying assumption we make is to take the diffusion constant of left- 
handed and right-handed quarks to be equal. Although this is a very good approximation 
since the diffusion properties are dominated by the strong interactions, it is one which must 
be treated with caution in certain limits of extreme suppression of baryon production by 
the decay processes. This is discussed in section 6 of [6] and the treatment given there can 
be applied to the present case. We will not discuss this here. 

With these assumptions we can now greatly simplify the full set of coupled equations 
(43) for all species in the plasma. Defining the variables 

fit = N c (jj tL - jj tR ) jj A = N c Hi(fi\ L - ji\ R ) jj B = N c E t (fil L + J (45) 

where ftt L means the difference in the chemical potential of the left-handed top and its 
(right-handed) anti-particle etc., the sum is over flavors and N c is the number of colors, 
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we can extract a simple set of three coupled equations: 



d Bl + Vw h _ r El _ Ef£^ = --f\z) - 2N c r f v w9A - 

D& + v w ^- -T ss ^ = -^F'(z) - 2N c T fVw9A ^ (46) 

J-o o -L o o 1 o 

T W T Q S T S T 
where F = 4cN c v w gA - ■ Nt is the number of fermion flavors - and we see that the 



first two equations reduce to a single one in the case Nf = 1. Ff, T ss and T s are the 
rates for the helicity-flipping, strong anomalous processes and weak anomalous processes 
respectively. With this new convention equations (46) take the form of diffusion equations 
corresponding to those in the companion paper [6], and the rates V are identical to the 
rates used there. T ss and T s can be read off from [6], and the helicity flip rate is derived 
in Appendix B: 



1 m? 2 rn 1 m t 



2 T 2 s 100 T 



r ss =^K ss aiT*K 8S ^ (47) 
T s =9N FKs aiT w k 



3 x 10 4 



where ^K ss a A s T A and K s a 4 T 4 are the number of strong and weak (respectively) anomalous 
processes per unit volume per unit time A^r is the number of fermion families. 

Before proceeding to analyze these equations and calculate the baryon asymmetry, we 
stop and review the numerous assumptions we have made in deriving the equations (46). 

• Assumption 1: L > ^ so that most particles in the plasma are indeed accurately 
described by the WKB approximation. Typical wall thicknesses are L ~ 20/T. 

• Assumption 2:v w < ^fj, the 'thermalization' condition for the applicability of our 
fluid ansatz. We calculated D 5/T, so for typical thick wall L ~ 20/T this is an 
extremely good for any wall velocity. 

• Assumption 3: v w < A=. We work to linear order in the wall velocity assumed 
smaller than the speed of sound in the plasma, so that perturbations can propagate into 
the region in front of the wall and give 'non-local' baryogenesis. This restricts us to modest 
wall velocities. 

• Assumption 4: T s < v^/D, so that the back-reaction of the baryon number violat- 
ing processes on the perturbed quantities can be neglected. Using (25) we see that this 
corresponds to v w > lO -2 ^^. Numerical simulations indicate k s ~ 0.1 — 1 so this is very 
consistent with favored wall velocities v w ~ 0.1 — 1. If T s > v^/D the sphaleron rate can 
equilibrate in front of the wall. 

• Assumption 5: All decay processes involving Higgs particles on the external legs can 
be neglected. This can be revised along the lines discussed in [6], and should lead only to 
minor numerical corrections. 
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• Assumption 6: The diffusion constants of quarks of opposite chirality can be taken 
to be equal. A very good approximation which need only be revised for the limit of extreme 
suppression of an asymmetry by decay processes. 

• Assumption 7: Only the top-quark Yukawa coupling is non-zero, so we do not 
describe the case of a two doublet model in which the lepton Yukawa coupling may be 
large. 

6. The Baryon Asymmetry 

The baryon number density B = jj-^UB, and calculating with the last equation in 
(46) we find its value at the front of the wall to be 

i i r r°° 



This is the final value of the baryon number in the broken phase under the assumption 
that the baryon number violating processes are turned off everywhere in the broken phase. 

2 

The expression (48) is valid with the assumption we made that T s satisfies » 1. 

To determine the final baryon asymmetry we must use the first two equations in (46) 
to extract /2a i n front of the wall which is the effective source for the baryon production. 

We will not solve exhaustively the two coupled equations for p, t and ft a m (46), but 
limit ourselves to the case that the dominant decay process is T ss so that we can drop the 
term Tffi t in the second equation in (46). This is justified for all k ss in the range 0.1 — 1, 
since < T 2 /4 on the wall (cf. (47)). We are then left with a single equation for fiA- 

D^ + v w ^-T^ = -^F'(z) - 2iV c r> w « M | (49) 

To solve (49) we need to find particular solutions as well as solutions to the homoge- 
neous equation. The latter are ~ e~ Xz where 

\d = ^(1 + ^ + ..) and -Ad = -£Kl-^ + ..) ^«1 

\={ r> w - (50) 

±A± = ±[(V% + -)TfB-] ^£ >>l 

The behaviour of the solutions is determined by the parameter — ^- (precisely as in the 

companion paper [6]). It characterizes the competition between decay and diffusion, 

being the time a typical particle spends diffusing in front of the wall before being caught, 
z^r 

When — ^ « 1, decay becomes irrelevant in front of the wall, and the only effect of 

V w 

the decay processes is to restore thermal equilibrium far behind the wall. This is also 
the criterion we need to use to see if any process is of relevance to the problem we are 
considering - in particular we used it above in (44), and implicitly in assuming that the 
Yukawa couplings which we have set to zero mediate decay processes which are slow in 
precisely this sense. 

We now consider the solution of (49) for each of the two source terms on the right 
hand side separately. 
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6.1 Classical Force Baryogenesis 

We return again to the 'ramp' ansatz (32), with A = 4cN c , and take only the force 
term as a source in (49). The problem is now homogeneous with boundary conditions 



D(-) +v w { — 

1 o J-o 



±- 



DF n 



o 



DLF n 



(51) 

derived by integrating (49) with the assumption that the integral of fx a is continuous across 
the boundaries. Requiring /I a to be finite at ±oo we take 
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x f z + B w e XbZ 



z > § 
z<-% 



(52) 



where Af, A w , B w and B b are constants, and from (50) we have 
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(53) 



Using the boundary conditions we determine the solution in front of the wall to be 
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(54) 



where X° D = ^ 



In various limits this reduces to 
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(55) 



where Z = Z{—^) and m D = m(— ^). 

The first two cases agree precisely with what we saw when we analysed the force 
neglecting decay processes in section 5.2. The prefactor in (55) corresponds to the average 
of the solution approached on the wall for << 1, but with the opposite sign so that 
the integrated contribution in front of the wall cancels that on the wall. The two limits 
for ^fy- are the limits of good and poor transport which we discussed. 

The two other cases given in (55) tell us how these solutions are modified when 

— 2^ > 1. The penetration of the diffusion solutions into the unbroken phase is reduced 
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since y > As we noted this is just the condition that the average diffusing particle's 
time in the unbroken phase before capture be longer than its decay lifetime. A second 
parameter enters in determining how the amplitude of the solution in front of the wall is 

/ r 

changed. In the first case, y -g-L « 1, the amplitude in fact compensates by increasing 
so that the integrated result jT p,A is unchanged (up to a factor of two). In the second 



case, 



^jL » 1, the amplitude is attenuated and the integrated result differs from 



that in the no-decay case by a factor of v 2 u /2DT^ « 1. The physical meaning of the 

parameter \J^g-L is also simple. It is (the square root of) the ratio of the time a particle 
takes to diffuse across the wall to its decay lifetime. So what the the third and fourth cases 
in (55) tell us is that the net density of particles in front of the wall in these stationary 
solution is not changed (up to a factor of two) unless the decay process is fast enough so 
that particles can decay as they cross the wall. This is a surprising result as one might 



expect the only relevant parameter to be 



DT, 



which compares the decay time to the time 



a particle actually spend in front of the wall. 

The factor of two has a simple explanation which will be familiar to the close reader 
of [6]. If one considers a diffusion/decay equation for a single species with a given injected 
flux modeled by a delta function, one finds two distinct regimes corresponding to the value 
of In the high velocity diffusion regime (^ << 1) the stationary solution puts most 
of the injected flux into the amplitude in front of the wall; in the low velocity case the 
solution has equal and opposite amplitude in front and behind, sharing the injected flux 
so that the amplitude is exactly half that in the first case. 

One might worry that some of these features are artefact of our ansatz (32) in which 
there is a discontinuity in the profile of m 2 Z at the back of the wall. To check this we take 
instead the 'bell' ansatz: 
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otherwise 



(56) 



and using the same conventions as in (52) we find in front of the wall 
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where \°, = — and A° n = ^. In the same limits as before this reduces to 
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(58) 
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This again shows the same results in each case. The only difference is the numerical factor 
which comes from averaging the profile over the wall, which is precisely what we anticipated 
since ?f =< Av w ^- >. 

Using (48) we now finally calculate the baryon asymmetry in its standard form and 

find 



n B 
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K S (X W 

9* 



9AZm? 
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dz 



if v 2 w < ^f - and \j ^ff-L ss 
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r s J? if v l < ¥ and 
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14 
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< 1 



(59) 



> 1 



where s = ^r-g*T 3 is the entropy density of the universe, g* the number of relativistic 
degrees of freedom, and L the thickness of the wall, rj is a geometrical factor which must 
be calculated for the particular wall profile (1 for 'ramp', 3 for 'bell'). In all cases we have 
assumed T s << v w /D; in the case of a wall moving sufficiently slowly that this condition 
is violated, one would recover ub oc v w as expected. Recall that these results are derived 
under the assumption that the strong sphaleron is the dominant decay process everywhere. 
Other cases can be treated using the more general form of (46) . 

This expression is remarkably simple. Most strikingly L, v w and D all cancel out in 
the answer in the most interesting (and quite plausible) regimes. For typical 'slow thick 
wall' parameter values e.g. L ~ v w ~ 0.1, k ss ~ 0.1 — 1 and rat ~ T, the conditions 
for our derivation hold. The result is « 2 x 10~ 2 k s o4(^) 2 A(9 2 x 10~ 8 k s A6>, where 
A9 = m~ t ~ 2 J m t (z) 2 g aZ dz is a measure of the CP violating condensate on the wall. 

The magnitude of A9 depends on the precise profile of the wall, with the greatest 
effect occurring (in the two doublet theory) if the phase 9 rolls fastest where the mass 
is large. In two doublet theories it will have the same sign on every bubble in a way 
determined by the effective potential, and can be 0(1) consistently with measurements of 
CP violation. In the case of the standard model Z condensate this factor will contain a 
suppression (potentially many orders of magnitude) depending on exactly how one sign of 
the condensate comes to dominate over the other. A full treatment of this case is required 
which goes beyond the scope of this paper [13]. 

6.2 Local Spontaneous Baryogenesis 

In the previous section we have concentrated primarily on the chemical potential 
fluctuations in front of the wall where the sphaleron rate is unsuppressed. This clearly 
dominates baryon production in the case of efficient transport ^fr < 1. In the case of 
inefficient transport (when the wall is very thick and slow) so that >> 1 the non- 
local baryogenesis will be suppressed by ~ and local baryogenesis may dominate. 
We analyse this case now to relate our treatment to that in the literature (in particular 
in CKN's work [8]) prior to [10], in which the potential importance of transport in the 
plasma was noted. 

When we combine the first two equations in (46) and ignore transport (D = 0), we 
obtain 

Vv/jlI -(T f + T SS )^' A + - —) f i A = -2N c T f v w g A Z' (60) 
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The force term source drops out in this limit simply because in order to induce perturba- 
tions particles must move in response to the force, (and D = "freezes" the particles). 
The case of the perturbations induced by the second source terms on the right of (46) is 
quite different. The effect of the field Z is local - it creates a perturbation at the point at 
which it is turned on by changing the interaction rates of hypercharge violating processes. 

For comparison with the previous literature it is instructive to consider the following 
cases in which the solutions to this equation can be read off simply. 

• Case 1: FJ 1 « L/v w « T" 1 for which 

A* A = 2N c v w g A Z (61) 

to leading order in ^j. 

• Case 2: T" 1 , FJ 1 « L/v w for which 

/*a = -2 ^ i fr-v w g A Z' (62) 

1 — iV j J- SS 

Here we ignore the homogeneous solutions which simply describe how the perturbations 
induced on the wall decay away behind the wall. 

The first case gives us what would be obtained by finding the local thermal equilibrium 
subject to the constraints imposed by the interactions locally, neglecting strong sphalerons. 
This is precisely the limit calculated by CKN in [8], albeit with a fermionic hypercharge 
in place of v w Z and correspondingly a fermionic hypercharge violating process in place of 
F f. The requirement FJ 1 << L/v w is just the condition that the interaction time for the 
hypercharge violating process be short in comparison to the time of passage of the wall, a 
requirement imposed by CKN on the fastest fermionic hypercharge violating process. 

From the second case we see that in if one takes T ss — > oo i.e. puts the strong 
sphalerons into equilibrium, the result is zero. This is a simple way of seeing the result 
obtained by Giudice and Shaposhnikov in [14]. The result (62) tells us how to calculate 
corrections to this constrained local equilibrium calculation, keeping the constraints but 
taking the rates to be finite, and gives what one might guess: the result is approximately 
equal to the equilibrium result of (61) with an additional suppression 

These formulae also show - as remarked earlier and originally pointed out by Dine 
and Thomas in [9] - that the perturbation is not well modeled by a fermionic hypercharge 
potential. But is it a potential for total hypercharge? In the rest frame of the plasma, with 
our assumption of a stationary wall profile, the time component of Z^ is v w Z. But one 
must be careful about the conclusion that one simply replaces v w Z by whatever this time 
component is. One evident problem is that one loses the explicit v w dependence, as such 
a potential can in principle be non-zero with the wall at rest, and the result when v w — > 
is then non-zero. The problem is exemplified in CKN's new calculation of spontaneous 
baryogenesis [12] which incorporates the effect of transport through a diffusion equation, 
in which only the time component of Z^ appears. However if we consider the case in which 
only the time component of Z^ is nonzero, their diffusion equation does not approach 
the correct thermal equilibrium as v w ^ 0. To recover this must include the force term 
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which results in a net overdensity in proportion to Z Q . In general, a correct treatment 
should include both space and time components of the gauge field, and the corresponding 
force terms. In our case, where the field configuration is assumed static in the wall frame, 
one can work in this frame and in this case the only non-zero component is the spatial 
^-component, and the time component plays no role. We will not compute in this case the 
baryon asymmetry, as it involves making specific assumptions about how the electroweak 
sphaleron rate behaves on the wall. We concentrate instead on the non-local variant of 
this mechanism. 

This brings us to one final remark. When the transition proceeds, as we have assumed, 
by bubble nucleation, is a spacelike vector if we take the wall profile to be stationary. 
The one case in which it can be modeled consistently as a timelike total hypercharge 
potential is when the transition occurs by spinodal decomposition, where the Higgs fields 
all roll together in the same way everywhere in space. Going back to our dispersion relation 
in section 2 for this case we would follow through our derivations in the same way. We would 
need to redefine the chemical potential by g^Zo and, because of the spatial homogeneity, 
would find no force term, and for the same reason, would discard the spatial gradients and, 
with certain assumption about how the transition proceeds, the time derivatives too. We 
would arrive at (61) above with v w Z replaced by Z Q . 
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6.3 Non-Local Spontaneous Baryogenesis 

We now turn to the case where this local effect is turned into a non-local one by the 
effects of transport. Perturbations can then be generated in front of the wall where the 
sphaleron rate is unsuppressed, in contrast to the local effect which has the problem that 
baryons must in the region where the effect which is vev squared suppressed is turned on. 
For a discussion of this effect we also point the reader to [12]. 

We turn again to (49), again assuming F ss to be the dominant decay process, taking 
the second source term and the 'ramp' ansatz: 

T,Z r _l <z< l 



To 1 otherwise 

The vev squared dependence in the rate is absorbed in the ansatz. With this ansatz /2a is 
then continuous everywhere but the equation is not homogeneous. 
The particular solution for this ansatz is 

^ = \ Vw r;^rrl— + T^L\ ~2< z <2 (64) 
1 ° I otherwise 

and the solution in front of the wall 

-^^(vFaoiK- 1 + 1 1 + (1 - 1 + ^'^'^ < 65 > 

using the same conventions again as in (52), and = • 
Again we simplify this in various limits to 



+2N c v w = — 



-^(^)e-*^-*) ^«1 and ^ « 1 
J> (££i)e-*<*-#) and ^«1 

v w L v vf„ ' D v± 



ly/%Le-V y/%L«l and » 1 

i fp; e -yfifi(*-it) a/%L»1 and ^»1 



(66) 



These solutions are very similar to those in the case of the force, the amplitude being of 
the form of the amplitude of a 'static' solution (61) on the wall multiplied by a factor 
which depends on the relative importance of decay and diffusion. The 'static' solution in 
this case is attained in the limit that there is no transport and the decay processes are 
turned on. 

In contrast, the solution (36) for the force is reached when transport is perfect and 
the decay processes are turned off. 

For the 'spontaneous' effect, we see that as D — > the amplitude in front of the wall 
vanishes rapidly (~ D 2 ). 

As D — > oo the amplitude goes as ^= but the tail integrates to give exactly half 

the solution (61), but in this case not with the opposite sign. The arguments which we 
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used to explain the sign in the case of the force do not apply as they rested on (40). In 
the 'spontaneous' result the sign is in all but the first case the same as that of the 'static' 
solution. What is happening is that the overdensity which the local process creates diffuses 
out to the region in front of the wall. Only in the first case, where the transport is much 
more efficient than the decay do particles diffuse in to 'cancel' the density on the wall. 

The baryon to entropy ratio can again be calculated. We find precisely the same 
expressions in (59) but with the replacement 

^ dz _ » f lll^ iz |^-^r < » r (67) 



7? 2c J r„ [l 

Using Ff and F ss from (47) so that Tf/r ss ~ 0A(mt/T) 2 / n ss , this result is seen to be 
equal to the force result multiplied by ~ 1 / 50 K ss a 2 ~ 1/k ss . (This is valid for « ^f- 
and, as we assumed, T ss » T/.) The final result is that the two effects have quite different 
parametric dependence but appear to be roughly of the same order of magnitude. 

A striking difference between the classical force effect and the spontaneous baryoge- 
nesis effect is in the dependence on the strong sphaleron rate. In the latter there is an 
inverse dependence on the sphaleron rate for a sufficiently large k ss . This dependence has 
been noted by CKN in their numerical study [12]. In contrast, the force sourced result is 
not suppressed by the strong sphalerons if ^k ss LT '/14 < 1. 

7. Comparison of Thin and Thick Wall Regimes 

It is interesting to compare the results in this paper with those in [6] , where we consid- 
ered the case of baryogenesis produced by reflection off a thin wall. The two calculations 
are most easily compared by looking at the version of the fluid equations (46). Once we 
established that the thermal fluctuations were unimportant we were able to reduce the 
system to a diffusion equation precisely analogous to that in [6] , but including a force term 
which extends over a finite region of space - the wall - where we had before the derivative 
of a delta function, modeling the injected flux. 

iJ Q 5 - jF(z) (68) 

In the classical force case, integrating the right hand side gives zero because the force is 
derived from a potential. Unlike the injected case, the effect of the classical force vanishes 
as the wall thickness is taken to zero. This is because the force simply speeds up and slows 
down particles that pass over it. We have not included the effect of the force as a reflecting 
barrier on low momentum states, since we neglected the effect of the 'flow' of one branch 
of the dispersion relation into the other at low momentum. If we include this we should 
obtain in the thin wall limit what one would calculate for the reflected flux in the WKB 
limit. The quantum mechanical reflection we can of course not recover. To compare the 
magnitude of the classical force effect with the WKB reflection effect we take the ratio of 
the amplitudes of the diffusion tails at the front of the wall. Doing this we find, neglecting 
decay and assuming efficient transport, 

J^IL „ ^-LT^-r, (69) 
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The most striking result is the the very different parametric dependence. 

The effect of introducing a decay term - in particular due to strong sphalerons - is 
quite different in the two cases. In our thin wall calculations [6] the suppression which 

2 

resulted was for r Vw n << 1, and this same suppression is seen in the nonlocal 

spontaneous effect. In the case of the classical force this suppression (which comes from 
the shortening of the diffusion tail in front of the wall) is compensated for by an increase in 
amplitude until T ss enters the regime \/T ss / DL > 1. An explanation may be found from 
the equation for /2a, (46). Dropping both Tf terms, integrating once one finds that the 

integral of /2a over all space must be zero. What the condition ^ ^-ff-L < 1 means is that 

strong sphalerons have little effect on the wall. But because J /j,a = 0, the compensating 
tail in front of the wall is also unaffected. No analogous conservation law holds in the 
spontaneous baryogenesis effect, because the whole effect is driven by a decay process (Tf) 
and strong sphaleron suppression occurs. 

8. Conclusions 

In this paper we have developed a new procedure to describe perturbations produced 
by a CP violating bubble wall moving through the plasma. 

We have shown how a Boltzmann equation can be used to describe the dynamics of 
thermal particles in the plasma when scattering processes are important. We now conclude 
with several remarks: 

• The treatment relies on a WKB approximation which is good for most of the particles 
in the plasma. It does not describe the effect of the background on low momentum particles 
for which this WKB approximation is not valid. A full treatment of all the perturbations 
produced by the propagating wall remains an open problem - for recent attempts see [15]. 

• We have neglected the effect of screening. The effects of the long range fields may 
be incorporated through the appropriate terms in the Boltzmann and fluid equations. We 
have not done this as it greatly complicates our analysis by coupling all species. However, 
in section 7 of [6] we discuss this issue in some detail. As illustrated by a simple model 
calculation of screening by quarks and leptons presented there, we do not expect that the 
effects of screening would alter our final results by a factor very different from unity. 

• There are many improvements and refinements of our calculation which are possible. 
The Boltzmann equation is in principle soluble without any truncation, and certainly there 
are other approximations which can be used. • One caveat must be added to our justifica- 
tion of the neglect of the Higgs particles in the determination of the perturbations driving 
baryogenesis. By dropping the Higgs particles we are assuming that they themselves are 
not significantly perturbed by a force term. In the thin wall case we also assumed that 
there was no injected flux in Higgs particles. This is not justified as the dynamics of the 
Higgs particles are likely to be non-trivial in the background of a changing vev. It is quite 
conceivable that such an effect could be important. 

• The methods we have developed should be useful in approaching the problem of 
the determination of the wall velocity and backreaction on the wall due to reflection. In 
particular we derived a set of force terms due to the mass which we did not make use of 
in the calculation reported here. 

• An important outstanding problem is to relate the techniques developed here to 
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some formal field theoretic methods such as those used in [16] and more recently in [17]. 
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Appendix A. Collision Integrals for Boson Exchange^ 

We consider first the contribution to the collision integrals in the fluid equations 
which come from the t-channel gauge boson exchange processes shown in Figure 2.° We 
will treat not just the quarks but also the right and left-handed leptons, primarily because 
it is interesting to compare our result for the diffusion constant with those obtained with 
our previous method in [6]. 

3 

When we integrate the Boltzmann equation as described in the text over J ct p, 
J p ct 3 p and J Pz^p there are six integrals we need to evaluate: 

/ " f P )fk(l ~ f k )\M\ 2 (2nf5\p + k-p' - k')(p -p' ) x \ Po 

Jp,k,p' ,k' 



Pz 

[ U(l ~ f P )fk(l - f k )\M\ 2 (2n)H\p + k-p>- k')( Pz - p' z ) x J Po 

Jp,k,p',k' { p z 



(70) 



where we have ignored the higher order terms in the perturbations ST, v and [i so that 
fermion population densities (in the plasma frame) are f Pi = 1/(1 + exp(p°/T )), pi = 
{p,k,p',k'}. 

For the gluon exchange process the scattering amplitude is 

l^l 2 = ^(^ ™ 

where A q = 32gj, t = (p — p' ) 2 ~ —2p-p', s = (p + k ) 2 ~ — 2p- k, u = (p — k' ' ) 2 ~ — 2p- k', 
g s is the strong coupling constant and m g is the thermal gluon mass. Note that here and 
in the following we neglect masses of the particles at legs of the diagram. Since typical 
scattering particles are thermal p ~ k ~ T, including them would introduce mass-squared 
corrections to our collision integrals which we can safely neglect. For the SU(2) and U(l) 
boson exchange diagrams for leptons 

lMwl2 = Aw jr^r |a1b|2= ^ b (F^if (72) 



t We thank Guy Moore for pointing out an error in an earlier version of this appendix. 
° In this paper we do not calculate all of t channel tree level diagrams that contribute to 
the diffusion constant. We refer the reader to [5] where this has been done systematically. 
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where Aw = 36g^, Ab = 78g^ tan 4 9 W , g w the weak coupling constant, 9 W the Weinberg 
angle, and mw and mg the thermal masses of the 577(2) and U(l) gauge bosons respec- 
tively. These amplitudes include the counting over all the fermions and anti-fermions which 
the given fermion can scatter off. Since the t-channel processes are dominated by the in- 
frared exchange momenta we have used the same approximation to the boson propagator 
as in Appendix C of [6] : we have assumed the exchange boson is Debye screened and hence 
the propagator is corrected by t — > t — m 2 , where m(T) is the corresponding Debye mass. 
This approximation is reasonable for longitudinal bosons; but for transverse bosons, since 
they are not Debye screened at the one-loop level, it underestimates the scattering rates. 

To see that the first and fourth integral are zero we simply use the symmetry of the 
matrix element under exchange of initial and final momenta i.e. p p' and k — > k'. 

Using parity invariance of the amplitudes {pi — > —pi) one can show that the third and 
fifth integrals are also zero. We are left with only two integrals which we denote as 7° (the 
second integral) and 7 3 (the sixth integral) which can be written in a Lorentz invariant 
form as follows: 



o r a s p r a 3 k o 3 r a s p r a 3 k 3 

1 =J 2p~J 2K fpfkL ' 1 =J 2p~J 2K fpfkL 

L " = J ~ fp ' ] I " ^l^l 2 ^) 4 ^ + k-p'-k>) 1 - [-t + {(p-p') ■ uY] 

(73) 

where is the plasma vector which is in the plasma frame w M = (1,0), and f Pi = 
1/(1 + ex.p(pi ■ u/T)) are the population densities. Dropping the Pauli blocking factors 
for the out-going states we can evaluate the integrals over p' and k' in the center of mass 
frame. In this frame = (7,77*), which can be written in terms of the plasma frame 
quantities as 7 = (pg z + k^/VsP 1 , = (j? 1 + k pl )/Vs^, s pl = 2p pl ■ k pl . For convenience 
we choose the axes such that p is along the z-axis: p = p(0,0, 1), v has zero azimuthal 
angle: v = v(sm/3, 0, cos/?) and p 1 = p'(sm9' cos(p\ sin9' sin^', cos^'). After some algebra 
and keeping only the leading logarithm, which occurs for small 9', we get 



A i ■ m2, V 

A Av 2 
3t 3 = — \hvps\ivl3) 2 + 2p 2 ] ln-V 
8tv l j mj 



(74) 



where A = 2A q for quarks and A = Ayy, A = Ab for the leptons. In order to recast this 
in a Lorentz invariant form, note that the Lorentz scalars we can construct from w M , p^ 
and W are u -p, u ■ k and p ■ k (since p-p~0, k • k ~ 0, u ■ u = 1). The unique solution is 
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(-fvpsmfl) 2 = u ■ p u ■ k — ^p, so that we have 
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(75) 



It is now possible to do the remaining integrations using the approximation J f x x n In xdx 
ln(n + 1) / f x x n dx, f x = (l + exp we find 



I u = I 6 = 



1024tt 5 



9Cf ^ 



rr?/ 



Now going back to the definitions of Ft and F v in (21) and (19) we have 



(76) 



rv = 



4/9 T ' 



r 



r„ = 



(77) 



where p a 



21 

8tt : 



C 4 Tq . For the quarks this gives, when we use m? = 8na s T 2 ss 3.6T 2 as 



described in Appendix A of [6] for a s = j at T a ~ lOOGeV, 



18a ?1 9 

r„ = 3r T = ^c* 2 m- — 



7( 4 n s 2yra s w 20 
Using the relation D = b/3aT v from (22), this gives us the quark diffusion constant 
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(79) 



For the leptons we calculate 
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jva w tan 2 6*^/ 



290 



(80) 



where we used m| = fa,„ tan 2 ^vk^ 2 ~ 0.04T 2 and m 2 w = ^a w tan 2 T 2 ~ 0.7T 2 , and 
taking a w = ^ and sin 2 6>vk = 0.23. 

These values for the diffusion constants agree very well with the values we obtained 
in Appendix A of [6]. The methods employed differ in that the first allowed a more 
general perturbation to the phase space density, but required an assumption about the near 
equality of the energy of ingoing and outgoing scattered particles. The method presented 
here is based on a more restrictive form for the distribution function, but involves no 
additional approximations. 

Appendix B. Collision Integrals for Decay Processes 
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In this appendix we evaluate the gluon exchange helicity-flip rate which dominates 
the hypercharge violating interactions. The relevant Feynman diagram is shown in Figure 
2. The top/Higgs helicity-flip processes (hypercharge violating Higgs exhange and Higgs 
absorption or emission) are slower because there are fewer particles to scatter off. 

The rate for a top quark of helicity A to scatter into one of helicity A' is 



12 

A^A' -yg 



f^fo(Po/T)^f (k /T)I f 

J 2 Po 2k a (gl) 



Wo 2 K 

on the wall where particles acquire a mass. The factor 12/T 3 = l/3noa is in accord with 
the rate definition r / in the diffusion equations (46) . Even though the masses are spatially 
dependent, for our purposes it is sufficient to evaluate the rate assuming a constant mass. 
The corrections to this approximation are of order IVm where / ~ 1/M g is the Debye 
screening length of the gluon so that l\Vm\/m ~ l/g s TL « 1. In (81) we consider an in- 
going fermion with mass mi momentum and helicity {p , A} which scatters off a fermion m^-, 
k (with arbitrary helicity) into {p' , A'}, k' with the appropriate spin dependent scattering 
amplitude .Ma^a'- The notation we use in (81) is familiar to the reader: no = 3(sT 3 /47T 2 , 
= 1.202, f (x) = l/(expx + 1), ct 3 k = d 3 k/(2n) 3 . We ignore Pauli blocking factors. 
In order to evaluate |.Ma^a'| 2 we define, following e.g. [18], the spin-four vector as 

where A is the helicity, so that s ■ p = and s ■ s = — 1. The helicity projection operators 
for a massive spin-1/2 particle eigenspinor u{p, A) are 

u(p, \)u(p, A) = + 75^ + m) (83) 

and similarly for anti-particles. 

The scattering amplitude squared reads 

IMa^vI 2 =- ^ ■ T 2 

[(p-p0 2 -m|] 2 16 

=Tr [(# + m)Y(¥+m)^] ( 84 ) 
T 2 =Tr [(jt> + m) lv X - (l + ^ (p' , A')) (/+ m) lt }- (l + ^ (p, A))] 

where A s = 32g| = bl2n 2 al (a s = g 2 /4:iv = 1/7) includes the 12 quarks and anti-quarks 
that a quark can scatter off via the gluon exchange. Rather than evaluating |.Ma^a'| 2 fc> r 
each pair {A , A'} we can define the helicity-flip amplitude 



2 A-, 



(85) 
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and the no-flip amplitude 



\M. 



+ 1 = |^Wa— A'| (86) 



2 A -A 



which simplifies the problem considerably. We still have to plough through rather lengthy 
algebra to arrive at an expression for 7\ • T 2 ± 

1 11 

—T x -T 2 ±=- V — Ti • T 2 = 
16 2 f-J 16 

(p ■ kp' ■ k' + p ■ k' p' ■ k — m\ k ■ k! — m\ p ■ p' + 2m\m\) 

±- ^2 [~s ■ s' (p ■ kp ■ k' + p ■ k' p ■ k - p ■ p k ■ k') 

2 A'=A 



(87) 



— p • p' (k • s k' • s' + k • s' k! • s) + p • kp' • s k' • s' 
+ p ■ k' p' ■ s k ■ s' + p' ■ k p ■ s' k' • s 
+ p ■ k' k • s p • s' — k ■ k' p ■ s' p ■ s 
+ ml (k ■ s k! • s + k ■ s r k! • s — m\ s ■ s')] 

where s M = s M (p, A) and s'^ = s^{p', A'). 

Since this form is Lorentz invariant, the integration over the out-going momenta p' 
and k' is rather straightforward in the center of mass frame: p + k = = ff + k'. To 
leading order in m\ we obtain 



J_=^ln 
87T p ■ k 



1 + 



2p-k 

9 



A s p- k 



(88) 



with the gluon thermal mass squared Mg = 8jva s T 2 . 

We can now integrate both X_ (in the leading logarithm approximation as explained 
in appendix B of [6]) and X + according to (81) to obtain the helicity-flip rate T_ and the 
no-flip rate 

24 In 2 2 m? 2 



Tf =r_ = —^aiiT 

r + =— a * T 



n 3 T 2S^ 
*2 



(89) 



where i = J Q 2 ln(l + x\z/2Tia s ) ~ C2 + (1/2) ln 2 (xf /-Ka s ) — Txa s jx\ ~ 2.5 is the angular 
integral (z = p-k/pk) which is approximately unity; x\ ~ 1.3 is the value of the momentum 
p/T at which the momentum integrals in (81) peak. 

We use this estimate of Tf in Section 6. Using the same method one could evaluate 
the helicity-flip W- and S-exchange processes relevant for the leptons. However since our 
main focus in this work is on top quark mediated baryogenesis, we shall not do so here. 
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Appendix C. Finite Temperature Dispersion Relation 

In this appendix we follow Appendix B of [6] to arrive at a finite temperature Dirac 
equation in momentum space. In the presence of a Z field condensate it is convenient to 
write the equation in terms of 2 x 2 chiral spinors in the plasma frame 



(E - g A Z + c R )-a-(P- g A Z) 



y R + m R ^ L = 



m L V R + 



(E + g A Z + c L )+a-(P + g A Z) 



= 



(90) 



where the notation is that of [6]. We will consider a planar wall moving in the positive 
^-direction with velocity v w with a pure gauge condensate Z^ = {Zq, 0, 0, Z z ){z) so that 
Zq = —v w Z z . The dispersion relation we require is obtained by setting the determinant of 
(90) to zero: 



(E + cf 



P 2 -m% 



1 



1 



9aZ + - Ac ) + g A Z 



2 v -' 



4P'[g A Z + i-Ac) + 4g 2 A Z 2 [E + c) z — P± I — 8P z g A Z z (E + c) I g A Zo + ^Ac 



where we have defined 

cr + c L 



c = 



Ac = c L - c R , 



m 



ran 



(l + a L )(l + a R ) 



The dispersion relation for anti-particles is obtained by the replacement P z 
Zq — > — Zq (recall that Z z is even under CP). 

In the limit when both Zq, Z z — > we find for particles 



(91) 



(92) 



—P z and 



(E+ C y 



i 

— i 

2 



P\ ± ^Ac ) + m 2 T 



(93) 



which reduces in the high momentum limit to E 2 = P 2 + m\ + 2M 2 R . 

A second case which is easy to solve is when Z z = and Zq ^ 0. This is relevant to 
the case of spinodal decomposition (second order phase transition): 



P| t ( g A Z + ^Ac 



(E + c) 2 - 

which in the high momentum limit simplifies to 



-i 2 



+ m T 



(94) 



E 2 = 



\P\T9aZq 



m 



(95) 



Therefore, as long as P » Ml,r, the motion of a thermal excitation is not affected (to 
leading order) by the dispersive plasma effects. In a first order phase transition this is not 
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the relevant case since Z^ is space like; in the plasma frame for example \Z /Z Z \ = v w so 
for slow walls it makes sense to neglect Zq. 

To make a detailed comparison with the free particle case we now proceed to a more 
systematic study of (91) and rewrite it as follows 



A = E + c, 
= A 4 - 2 



P 2 + m\ + g 2 A Z 2 z + (g A Z + ^Ac^j A 2 + 8P z g A Z z [g A Z + ^Ac^j A+ 



P 2 + m 2 T - g\Z 2 z + g A Z + -Ac 



+ 4m| 



g A Z + -Ac 



(96) 

We now consider a weak Z field expansion of this equation and write the solution in the 
form: 

A = A (l + e) 



(97) 



where Ao is the solution of (96) with linear term in A neglected: 
A 2 = e 2 ± VAo 



e 2 = P z + m 2 T + g\Z z z + [ g A Z + -Ac 



1 



(98) 



A = 4 ( P 2 + g\Z 2 z ) ( g A Z + -Ac ) + 4g 2 A Z 2 (m 2 T + P 2 ) 



The leading order correction e reads 



2g A Z z P z (g A Z + \Ac) 



AoVAo 



(99) 



where the signs coincide with (98). In order to make more transparent what the dispersion 
relation (97) - (99) mean we look at its high momentum limit. 

We now restrict ourselves to the case Zq = 0, Z z ^ 0, which describes the case of a 
wall at rest in the plasma. Since for a stationary wall profile Zq = —v w Z z , we anticipate 
that this will be a good approximation for a sufficiently slow wall. There are two cases to 
consider depending on which term in the determinant A in (98) dominates 



VA ~2g A Z z ^p2 + m 2 T , for (Ac) 2 P 2 « 4g 2 A Z 2 (P 2 + m 2 ) 



E 



v / A ^|Ac|P~(M 2 -M 2 )-, for (Ac) 2 P 2 » 4g 2 A Z 2 (P 2 + m 2 ) 



(100) 



In the high momentum limit these two cases reduce to a comparison between Ml — 
and 2g A ZP z . Since M| - M\ a w T 2 (see [6] ), and we can write 2g A Z = Qcp/L, the 
two cases become approximately Qcp > Ci w LT and Qcp < oi w LT, which we therefore 
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refer to as cases of 'strong' and 'weak' condensates respectively. (Note that for a realistic 
wall thickness in order to be in the 'strong' regime Ocp must be of order unity.) 
In the first case the dispersion relations are 

(E L > R + cf = P 2 + (sign yjp* + m 2 T ± g A Z z ^j + (Ac/2) 2 

w) 7& 2(£l ' fl + c) + (Ac/2)2 ^ 

(e l > r + c) 2 = Pi + (sign P z yj P 2 + m 2 T Tg A Z^j + (Ac/2) 2 

T(Ac/2) vftr 2 ^" +c ) +(Ac/2)2 ^ 



(101) 



p2 
± z 

+ m% 



We have adopted here the notation used in the main text - labelling the WKB states 
by the chiral states they deform into in the symmetric phase. A careful inspection of 
these relationships in the high momentum limit P >> M LiR , reveals that they reduce, to 
leading order in coupling constants, to the zero temperature dispersion relations plus the 
finite temperature mass \[2M^ for the left handed particles and their anti-particles and 
\[2M R for the right handed particles and their anti-particles. 

Firstly we note that in the free case (when we set T = 0) (101) reproduces precisely 
the dispersion relations we had in (4) and (5). Secondly, when Z z = we see that the 
effect of the thermal corrections is to split the left-handed and right-handed states (by an 
amount oc Ac), but to leave the degeneracy of particles and anti-particles intact. Finally, 
and most importantly, we see that the corrections to the free dispersion relation are small 
in this case Ocp > Ci w LT, so that the analysis in the main text applies. 

The second case of a 'weak' condensate in (100) gives dispersion relations 



(E L > R + c) 2 = (P T Ac/2) 2 +m%± 2g A Z z ^ yV T Ac/2) 2 + m\ 



(102) 



(P z <« + c) = (P ± Ac/2) 2 +m 2 T T 2g A Z z ^yJ(P ± Ac/2) 2 + 
It is instructive to re-write the first relation in the limit when P >> Ml^r 

(E L > R ) 2 ~ P 2 + m 2 + 2M\ R ± 2g A Z z j r y / P 2 + m 2 (103) 

When m = there should be no physical effect due to the Z field since it is then just 
pure gauge. We see that this is the case since the dispersion relation can then be written 
(E L > R + c) 2 = (P'=f Ac/2) 2 where P' = |P'| and P' = (P ± ,P z ±g A Z z ) (compare this form 
with (94)). When rn^O however, just as in the zero temperature case, the dispersion re- 
lation (102) leads to a non-zero acceleration: v z = Td z (g A Z z m 2 )(P^ / P 3 E 3 ) - d z m 2 /2E 2 , 
where we have assumed that E and P± are conserved in the plasma frame, an approxi- 
mation correct to leading order in v w . Comparing this result with (8) we see that, in this 
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second case of a 'weak' condensate, the force term has the same form (oc d z (gAZ z m 2 )) but 
a somewhat different momentum dependence. Corresponding corrections to the analysis 
presented in the text would be required to describe this case, which we anticipate would 
lead to minor numerical changes in the coefficients of the terms in the fluid equations if 
the analysis is carried through in the same way. 
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Figure Captions 

Figure 1: Dispersion relation for particles in a background axially coupled pure gauge field. 

The plots show the energy E as a function of canonical momentum p z for p± = 0, 
for (i) m > ', (ii) m < gAZ and (hi) m = 0. 

Figure 2: Vector boson exchange diagrams which dominate in damping temperature and 
velocity perturbations in the fluid, and play a leading role in determining the 
diffusion properties of different particle species. The gluon exchange diagram 
also, in the presence of a quark mass term, describes the key hypercharge-violating 
process, namely the helicity flip process computed in Appendix B. 

Figure 3: Chirality flipping Higgs process which contributes to the damping of a chiral 
chemical potential in front of the bubble wall. 

Figure 4: Solutions for the chiral chemical potential fi in the background of a bubble wall, 
with a 'ramp' ansatz for the CP violating condensate field m 2 Z. The solutions 
for ^- are sketched for the two cases (i) « 1 and (ii) » 1. 
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This figure "figl-2.png" is available in "png" format from: 
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